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1. Introduction 

Copula is with no doubt the most popular tool for describing the dependence 
between two random variables. The popularity is partially based on the fact 
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that the dependence between any random variables can be modelled by some 
copula. This fact is known as Sklar’s theorem, which states that for any two 
random variables Y) and I 2 there exists a copula C (a two-dimensional real- 
valued distribution function with domain [0, l] 2 and uniform margins) such that 

P{W < u u Y 2 < u 2 } = c(p{Yi < m},P{F 2 < u 2 }), (1) 

for any u\,u 2 > 0. We refer to Cherubini et al. (2004), Joe (1997), Nelsen(2006) 
for a comprehensive overview of the copula theory. 

Now let us switch from random variables to stochastic processes and try to 
describe dependence between components of some two-dimensional Levy process 
X(t) = (Xi(t), X 2 (t)), that is, of some cadlag process with independent and 
stationary increments. Applying Sklar’s theorem for any fixed time moment t, 
we get that the dependence between X\ and X 2 can be described by some copula 
C t . i.e., 

P{Ar(t) < u u X 2 (t) < u 2 } = C t (p{Xr(t) < «i} ,P{X 2 (t) < « 2 }), (2) 

for any u\,u 2 > 0. Nevertheless, the direct application of the representation (2) 
to stochastic modeling has a couple of drawbacks. First, it turns out that the 
copula C t in most cases essentially depends on t, see Tankov (2004) for examples. 
Second, since the distribution of X(t) is infinitely divisible, (2) is possible only 
for some subclasses of copulas. In other words, the class C t depends on the class 
of marginal laws Xi(t). 

To avoid such difficulties, researches are trying to characterize the dependence 
between the components of Levy process in the time-independent fashion. One 
of the most popular approaches for this characterization is the so-called Levy 
copula (defined below), which was introduced by Tankov (2003), and later stud¬ 
ied by Barndorff-Nielsen and Lindner (2004), Cont and Tankov (2004), Kallsen 
and Tankov (2006), and others. Among many papers in this field, we would like 
to emphasize some articles about statistical inference for Levy copulas (mainly 
with applications to insurance), which include some basic ideas that are widely 
used in statistical research on this topic, in particular, in the statistical analy¬ 
sis in the current paper - Esmaeili and Kliippelberg (2010 and 2011), Avanzi, 
Cassar and Wong (2011), Bucher and Vetter (2013). 

The main objective of this article is the application of the Levy copula 
approach to a class of stochastic processes, known as time-changed Levy pro¬ 
cesses. In the one-dimensional case, the time-changed Levy process is defined as 
Y s = L'f(s), where L is a Levy process and T is a non-negative, non-decreasing 
stochastic process with T(0) = 0 referred as stochastic time change or simply 
stochastic clock. If the process T is also a Levy process, then it is called the 
subordinator , and the process Y s is usually refered as the subordinated process. 
The economical interpretation of the time change is based on the idea that the 
“business” time T(s) may run faster than the physical time in some periods, for 
instance, when the amount of transactions is high, see Clark (1973), Ane and 
Geman (2000), Veraart and Winkel (2010). 
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In this paper, we consider one natural generalization of the aforementioned 
model to the two-dimensional case. Our construction is based on the so-called 
multivariate subordination , introduced by Barndorff-Nielsen, Pederson and Sato 
(2001). In particular, we prove the series representation of the processes from 
our class, which allows to simulate the processes with given characteristics, and 
show an application of this representation to real data. 

The paper is organized as follows. In the next two sections, we shortly explain 
the notion of Levy copula and the idea of stochastic change of time. Afterwards, 
in Section 4, we introduce our model and discuss some properties of it. Our main 
results are given in Section 5, where we also provide some examples. In the last 
section, we apply our model to some stock prices. 

2. Levy copulas 

The construction of the Levy copula is based on the concept of tail integrals. 

Definition 2.1. For a one-dimensional Levy process Z with Levy measure Vz , 
tail integral is defined as 


U z (x) '■= 


vz ( x , +oo) , if x > 0, 
—vz{— oo, x), if x < 0. 

Definition 2.1 can be equivalently written as 


U z (x) := (-1)-^ (/(*)), 


where 


1 0*0 : = 


(x,+oo), 

(-00, x), 


if x > 0, 
if x < 0, 


and s(x) 


2, if x > 0, 
1, if x < 0. 


The reason for this definition is that in the case of infinite measure vz , Uz(A) 
is infinite for any set A which contains 0. Analogously, for a two-dimensional 
process Z = (Z±, Z 2 ) with Levy measure vg, 


Ug[x u x 2 ) := (-1) S(X1 > +S(X2) ^(/(*!) x I(x 2 j), 


and this definition is also correct for any real X\ and X 2 ■ 

— 2 — 

Definition 2.2. A two-dimensional Levy copula is a function from M to M 
such that 

1. F is grounded, that is, F (m, u 2 ) = 0 if m = 0 for at least one * = 1,2. 

2. F is 2-increasing. 

3. F has uniform margins, that is, F(u,o o) = F(oo,u) = u. 

4 . F (*ii, U2) 7^ 00 for (iq, u 2 ) ^ (00, 00). 
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The main result on Levy copulas is an analogue of the Sklar theorem for 
ordinary copulas which states that for any two-dimensional Levy process Z 
with tail integral U? and marginal tail integrals Uz 1 and Uz 2 , there exists a 
Levy copula F such that 

Ug{xi,x 2 ) = F (Uz 1 {x 1 ),Uz 2 (x 2 )) , (3) 

and vise versa, for any Levy copula F and any one-dimensional Levy process 
with tail integrals Uz 1 and Uz 2 there exists a two-dimensional Levy process with 
tail integral Ug given by (3) with marginal tail integrals Uz 1 and Uz 2 ■ The first 
part of this theorem can be easily verified for the case when the one-dimensional 
Levy measures are infinite and have no atoms, because in this case the Levy 
copula is equal to 

F( Ul ,u 2 ) = U{U^( U 1 ),U^(U 2 )), ( 4 ) 

where U is the tail integral of the Levy measure of X(t), see [20]. 

3. Time-changed Levy models 

As it was already mentioned in the introduction, the time-changed Levy process 
in the one-dimensional case is defined as 

Y s = L ns) , (5) 

where L is a Levy process, and T(s) - a non-negative, non-decreasing stochastic 
process with T(0) = 0. This class of models has strong mathematical background 
based on the so-called Monroe theorem [22] , which stands that any semimartin- 
gale can be represented as a time-changed Brownian motion (that is, in the 
form (5) with L equal to the Brownian motion W) and vise versa, any time- 
changed Brownian motion is a semimartingale. Various aspects of this theory 
are discussed in [5] and [11]. Nevertheless, the first part of the Monroe theorem 
doesn’t hold if one introduces any of the following additional assumptions: 

1. Processes W and T are independent. This assumption is widely used in 
the statistical literature and is quite convenient for both theoretical and 
practical purposes. 

2. Time change process T is itself a Levy process; such processes are known as 
subordinators. In this case any resulting process Y s is also a Levy process, 
which is usually called a subordinated process. 

These drawbacks of the time-changed Brownian motion lead to the idea of con¬ 
sidering more general model (5) with any Levy process instead of the Brownian 
motion and introducing the assumption that the processes T and L are inde¬ 
pendent. This model has been attracting attention of many researches, see, e.g., 
[7], [8] [10], [12], [26], _ 

Nevertheless, there is no clear understanding in the literature how to extend 
this model to the two-dimensional case. The most popular construction is to 
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consider the model (5) with a two-dimensional Levy process L and to provide 
a time change in each component with the same process 7”, see [25]. 

Interestingly enough, in the case when L is a Brownian motion, the correlation 
coefficient between subordinated processes is upper bounded by the correlation 
coefficient between the components of the Brownian motion, see [16]. Moreover, 
these coefficients coincides in some cases, see [15]. 

4. Two-dimensional subordinated processes 

In this section, we introduce a two-dimensional generalization of the model (5). 
This generalization is based on the notion of the two-dimenstional subordinator, 
which we define below. 

Definition 4.1. A two-dimensional subordinator T(s) = (7i(s), T 2 (s)) is a 
Levy process in 1R 2 such that both components T\ and T 2 are one-dimensional 
subordinators. 

Evidently, T is a two-dimensional subordinator if its margins Ti and T 2 are 
independent. Another example is given by the following statement (see [27]): if 
T\ , T 2 , T 3 are 3 independent subordinators, then the processes 

(Ti00 + r 3 (s), t 2 (s) + r 3 ( s )) 

and 

(Ti(T 3 ( S )), T 2 (T 3 (s))) 

are two-dimensional subordinators. To the best of our knowledge, general criteria 
providing necessary and sufficient conditions for two-dimensional process with 
marginal positive Levy processes to be a two-dimensional subordinator, are not 
known in the literature. 

Consider now a two-dimensional Levy process L(t) = (Li(t), L 2 (i)) with inde¬ 
pendent components and a two-dimensional subordinator T(s) = (7i(s), ^(s)) 
such that Ti{s) is independent of Z/j(s) for any i = 1,2. Define the subordinated 
process by composition 

X(s) = (x 1 (s),X 2 (s)) := (L^T [(*)),L 2 (T200)). (6) 

This construction, known as multivariate subordination, was firstly considered 
in [4], The next theorem sheds some light to the characteristics of such processes. 

Theorem 4.2. Let W l , i = 1,2 be two independent one-dimensional Brownian 
motions and let T(s) = (71 (s), 72 (s)) be a two-dimensional subordinator with 
Levy triplet (p,0,rj), where p = (pi,p 2 ) with pi > 0, i = 1,2 and p is a Levy 
measure in 1R 2 + . 

Denote by diag(a:, 2 /) with x,y € 1R a two-dimensional diagonal matrix with 
the values x and y on the diagonal. 

Then 
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• the process 

X(s) := (wMs^WMs))) (7) 

is a two-dimensional Levy process; 

• the Levy triplet of the process X is given by 

(o,diag(pi,p 2 ) 

where the Levy measure v is defined as 

v{B) ■= J 2 p[b\ diag(3/i,2/ 2 )) r]{dy 1 ,dy 2 ), B C B 2 , 

and pt(B]A) := P{£a £ B} is the probability that a random variable £4 
with zero mean and covariance matrix A belongs to the set B. 

Proof. This theorem is essentially proven in [4]. O 


5. Series representation for subordinated processes 


In this section, we apply the result by Rosinsky [24] to our setup. 

Theorem 5.1. Let X(s) be a two-dimensional Levy process constructed by mul¬ 
tivariate subordination of the Brownian motion, see Theorem f.2 for notation. 
Denote by F(u,v) a positive Levy copula between 7i(s) and T 2 (s). Assume 
that F(u,v ) is continuous and the mixed derivative d 2 F{u,v)/dudv exists in 
Bff. Moreover, assume that there exists a density function p *(•) and a function 
f*(u,x ) : B\ —> B+, such that 

1. for any u, x > 0 , 


r f*(u,x) 


p*(z)dz = 


dF(u, x) 
du 


( 8 ) 


2. the function f* (u, x) monotonically increases in x for any fixed u, and 
moreover the equation 


f*(u,x ) = y 


has a solution in closed form with respect to x for any y > 0 ; we denote 
this solution by h*{u,y). 

Next, define a two-dimensional stochastic process Z{s) = {Z^s), Z 2 {s)): 

OO - - 

Z ^ s ) ■■= Y,V U ■!{/*< a}, (9) 

i =1 

00 I - 

Z 2 (s) := Y, G i 2) J U t 1) (h*F i ,G[ 3) j)-I{R i <8} 

1=1 " 


( 10 ) 
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where U\ and U 2 are tail integrals of the subordinators T\ and T 2 resp., U[ ^ 
and U 2 1 '* are their generalized inverse functions, that is, 

U-~ 1 \y) = mf{x>0: Ut{x) < y} , i = 1,2, y £ M+, 

r, is an independent sequence of jump times of a standard Poisson process, 
Gp, Gp - are two sequences of i.i.d. standard normal r.v., GP - sequence of 
i.i.d. random variables with density p* (•), Ri -sequence of i.i.d. r.v., uniformly 
distributed on [0,1], and all sequences of r.v. are independent of each other. 
Then 

X(s) = Z(s), Vs €[0,1], 

Examples. 

1. Consider the positive Clayton-Levy copula 

Fc(u,v ) = (u~ e + v~ e )~ 1 * * ^ 6 

with some 9 > 0. First derivative of this function is equal to 

dFc{u,v) (y/u) e+1 

du ~ (l + ( v /u) e ) {1+9)/9 ' 

Motivated by this representation, we suggest to define the density function 
p*(z) as 

P * (Z) = dz {(i + ^)(w} = dz { (2 -* + l) (1+ *>/ e }’ 

and the function f* ( u , x) := x/u. Both conditions on the functions p* and 
f* are fulfilled. 

2. Note that the same arguments can be applied to any sufficiently smooth 
homogeneous Levy copula, that is, to any copula such that 

Fn(ku,kv) = kFn(u,v), \/u,v,k> 0, (11) 


see Remark 5.2 about the difference between ordinary copulas and Levy 
copulas. Note that taking the derivatives with respect to u from both parts 
of (11), yields 


d d d 

— F H (ku,kv ) = —F h (u,v) = —F H ( ri ,r 2 ) 


and therefore one can define 


>-1=1 

7*2 =v/lL 


p,(2) = s sV ? ' ,(ri ' r2) 


r i=l 
r 2 =z' 


f*(u,x) := x/u. 


For the description of the class of homogeneous Levy copulas we refer to 
Section 4 from [3]. 
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3. Moreover, we can apply the same approach for any mixtures of homoge¬ 
neous Levy copulas. In fact, consider the function 

n 

F m (u,v ) = y ^/3 r F r {u,v), 

r—1 

where /3i, ...,j3 n are positive numbers such that X^r=i fir = 1 and F r (u,v) 
are homogeneous Levy copulas for any r = 1 ..n. It’s easy to see that 
Fm(u,v ) is a homogeneous Levy copula. As it was shown in the previous 
example, we can take f*(u, x) := x/u. Note also that in the case of mixture 
model, 

n 

p*{z ) = '^2PrP*(z), 

r= 1 


where by p*(-) we denote the density functions such that 
rx/u dF r (u, x) 


/ X/U 

p*(z)dz = 

-OO 


du 


Proof, (of Theorem 5.1) Since the Levy copula F is sufficiently smooth, we can 
differentiate both parts in (3) and get that 


r , \ F 

<B> = L appr 


UK: 


p(b ■,Aitxg{U l 1 (ri),U. 


d(U 1 (y 1 ))d(U 2 (y 2 )) 

r-i=(7i(yi) 

r2=U2{y2 ) 

d 2 F( n ,r 2 ) J J 

—-— - dridr 2 , 

or i or 2 


see Proposition 5.8 from [15]. Our general aim is to apply the result by Rosinsky 
[24], which is nicely formulated as Theorem 6.2 in [15]. Comparison of this 
result and the statement of our theorem leads to the idea to find a function 
H : (0, Too) xS-> JR 2 , where S' is a measurable space, such that 


(. B )=J P< H (r,^) G B\dr, V B £ B(1R 2 ), 


( 12 ) 


where D is a random vector from S. 

First note that it is sufficient to consider the sets B = B\ x B 2 , where 
B\ = [x, oo), B 2 = [y, oo), x,y £ M. For such B , 

p(B\ diag(l/ 1 _1 (ri),t/ 2 _1 (r 2 ))) = ;U/ 1 {r 1 )^ ■ p(b 2 ; f7 2 _1 (r 2 )), 

where by p(o,B) we denote the one-dimensional normal distribution with zero 
mean and variance equal to o (since there is no risk of confusion, we use the 
same letter for 2-dimensional and 1-dimensional distributions). Therefore, 


v{B) 


J EjSf(n) p(b 2 ; U 2 1 (-)) 


dr i, (13) 








V.Panov and I.Sirotkin/Series representations for time-changed Levy models 


9 


where by 







dr 2 


we denote the mathematical expectation with respect to the measure Jz?(ri) with 
the distribution function F(r 2 ) = dF(ri, r 2 )/dri (the statement that F(r 2 ) is 
in fact a distribution function is proven in [15], Lemma 5.3). By the well-known 
Fubini theorem, 




(■)) 



g{v ;n) dv, 


where 


s(-; n) 



d 2 F(n,r 2 ) 

dridr 2 


dr 2 , 


n > 0, 


(14) 


and pi (• ; U 2 1 (?' 2 )) is the density of the normal distribution with zero-mean 
and variance equal to 1 (r- 2 ), and d 2 F(ri,r 2 )/(dridr 2 ) is the density function 
corresponding to the distribution function F{r 2 ). 

Note that g is a density function, see Remark 5.3. Changing the variables we 
get 


5(S r i) 



d 2 F(n,r 2 ) 

dridr 2 


d(U 2 (f 2 )), 

r2=U 2 {f2) 


r i > 0. 


The last expression yields that g(-', r\) is in fact a variance mixture of the 
normal distribution (see [6] or [21]). This in particularly gives that the random 
variable 


£ = V1VV2 


has a distribution with density g(-;r i), where rji has standard normal distribu¬ 
tion, and tj 2 - distribution with density 


P{",ri) 


d 2 F(r 1 ,r 2 ) 

dridr 2 


U&) 

r2=U 2 (-) 


d 2 F(r 1 ,r 2 ) 

dridr 2 


m-)\, 

r 2 =U 2 (-) 


that is, the density of the random variable U 2 where rjs has a distribution 

function F{r 2 ). Since (8) holds, we get that 


d 2 F(r 1 ,r 2 ) 

dridr 2 


df*(ri,r 2 ) 

dr 2 


P*{f*(ri,r 2 )), 
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and therefore r / 3 has the same distribution as h*(r 1 , 1 ) 4 ), where rj 4 has distribu¬ 
tion with density p*(-). 

Finally we get the following representation for the Levy measure w. 


"(B) 



U l 1 {r 1 )^jdu ■ P{ m yju 2 1 (h*(r 1 , 774)) 



dr 4 . 


This representation motivates to define the function H by 


H{r,D) 


f y Jup(r)-D 1 \ 
\D 2 y/u^(h*(r,D 3 ))J 


with D = (Z?i, D 2 , D 3 ), where Di,D 2 have standard normal distribution, and 
Z ?3 has a distribution with density function p*(-). This observation completes 
the proof. 


O 


Remark 5.2. In the context of ordinary copulas, it is common to introduce the 
homogeneous copula Cof order k by 

cP(ku,kv) = k a C H {u,v), \/k,u,v>0, (15) 

see, e.g., [23]. Substituting u = v = 1, we get Cn{k,k) = k a . Therefore, taking 
into account the Frechet bounds, we arrive at the inequality 

max (2A: — 1,0) < k a < k, 

which yields that a £ [1, 2]. Moreover, it turns out that the class of homogeneous 
ordinary copulas coincides with Cuadras-Auge family. More precisely, 

Cfl\u, v) = (min(it, v)) 2 a (uv) a 1 , u, »f[0,1], 

see Theorem 3-4-2 from [23]. Returning to Levy copulas, we realize that similar 
to (15) equality 

Fn(ku,kv) = k a FH(u,v), V k,u,v > 0 


is possible only in case a = 1 .In fact, taking limit as v —> 00 , we get the equality 
ku = k a u , \/u, which leads to trivial conclusion a = 1. This argument yields the 
definition of homogeneous Levy copula (11). 

Remark 5.3. Let us shortly show that the function g{r\] •) defined by (14) is 
a density function for any ri. In fact, as it was mentioned before, the function 
F{ r 2 ) = dF(r\, rfj/dri is a distribution function, and moreover d 2 F(r\,r 2 ) / dr\dr 2 
is the denisity function of this distribution. Therefore, g(r±; •) > 0, and 


/ g(rv, v)dv = / pi[U 2 1 (r 2 ); v) 

JiR JiR v 7 


dv 


d 2 F( ri ,r 2 ) 
dr id r 2 


dr 2 


[ d 2 F(n,r 2 ) 

/ r + dridr 2 


dr 2 = 1. 
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Remark 5.4. It is a worth mentioning that the right way to truncate series 
in (9)-(10) is to fix some r and keep N(r) = inf,; {T^ > r} terms, see [15] for 
details. 

6. Empirical analysis 

In this section, we consider the following model: 

X(s) := (^i(Ti(s)),IE 2 (T 2 ( S ))) (16) 

where 

Wi(t) = Hit + aiWi(t), * = 1,2, (17) 

Wi(i), W 2 (i) are two independent Brownian motions, /*i,/* 2 G 1R , eri,er 2 G M+, 
and (7i(s), 72 (s)) is a two-dimensional subordinator. The dependence between 
71 (s) and ^(s) is described via some Levy copula which belongs to a 

class parametrized by S G 1R. The marginal subordinators 7i(s) and 7 ~2(s) belong 
to some parametric class of Levy processes. The corresponding parameters are 
denoted by 9\ and 0 2 . 

In what follows, we will apply this model to the modeling of stock returns. In 
this context, X(s) represents the returns of two stocks traded on the Nasdaq, and 
(7i(s),7i(s)) are cumulative numbers of trades of these stocks. Our approach 
can be considered as a generalization of the paper [1] , where the one-dimensional 
time-changed Brownian motion is used for representing one-dimensional stock 
returns. In [1], it is shown that the cumulative number of trades is a good 
approximation of business time. More precisely, the authors showed that the 
theoretical moments of the subordinator almost perfectly coincide with the em¬ 
pirical moments of cumulative number of trades. 

Our simulation study consists of three steps. 

1. Estimation of the Levy copula between one-dimensional subordinators. 
First, we assume some parametric structure of the Levy copula between 
7i(s) and 70 (s) and estimate the parameters of this structure. Our esti¬ 
mation procedure is motivated by the research [17], [18] and is described 
below in Section 6.2. 

2. Estimation of the parameters of the processes W\(t) and IT 2 (t). On this 
stage, we apply methodology described in [1] to the processes W\{t) and 
W 2 (t), and get the estimators of the parameters /*i, /z 2 , <ti, ct 2 . This part 
of the empirical analysis is given in Section 6.3. 

3. Applying simulation techniques. Since we already have the estimates of the 
Levy copula between the subordinators and the parameters A*i, /u 2 , o - !, cr 2 , 
we can apply the main theoretical result of the current paper presented in 
Theorem 5.1. The algorithm is described in Section 6.4. Some graphs are 
given in the appendix. 
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6.1. Description of the data 

We examine 10- and 30-minutes Cisco, Intel and Microsoft prices traded on the 
Nasdaq over the period from the 25. August 2014 till the 21. November 2014. 
For each equity we have time variable, price variable and number of trades. 
The length of time series is 832 observations for 30- minutes data and 2496 
observations for 10- minutes data. Before the analysis, last observations of each 
trading day were deleted due to the abnormally small number of trades. We 
suppose that it is related to the microstructure of the market. 

Some descriptive statistics of the data are given in Tables 3 and 4 (see Ap¬ 
pendix A). 

6.2. Step 1. Levy copula estimation 

The techniques for the parametric estimation of Levy copula are not well- 
described in the literature and are mainly known for compound poisson pro¬ 
cesses (see [2], [17], [18]). In this research, we also assume that subordinators 
7i(t) and T 2 (t) are compound poisson processes (CPP) with positive jumps, 
that is, 


JVi(t) N 2 (t) 

7i(t) = Y, x it T 2 (t)=^Y j , (18) 

*=» 3 =1 

where X t and Yj are i.i.d random variables with densities fi(x; 9\) and f 2 (x', 9 2 ) 
with support in 1R + , Ni(t) and N 2 (t) are Poisson processes with intensities Ai 
and X 2 resp. 

At first glance, (18) seems to be a strong assumption. However, the Levy 
processes with truncated jumps, that is processes in the form 

jt = 'y ] A Z s , A Z s = Z s — Z s _, 

0<s<t 

\AZ s \>c 


where Z is a (multi-dimensional) Levy process, c is a positive constant, are com¬ 
pound poisson processes (see [15], [25]). Moreover, the truncation of the jumps 
is a well-used procedure for different simulation and estimation techniques, and 
is a quite natural tool in the context of stock data (see [15], [18]). 
Two-dimensional CPP could be represented in the following way: 

71(f) = T^{t) + rht), 

r 2 {t) = r 2 ± (t) + r^(t), 

where T^~(t), (7^(t), T 2 (f)) are independent compound Poisson processes, 
and moreover 
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• the Levy measure of the process T^~(t) is equal to ? 7 (-, 0 ), that is, the 
component T\~{t) represents the jumps of the process 7i (t) which occur 
independently of the process T 2 (i); 

• analogously, the Levy measure of the process is equal to 77 ( 0 , •), that 

is, the component represents the jumps of the process T 2 (t) which 

occur independently of the process 7j(t); 

• the Levy measure of the joint process ^7^ (t), T 2 (t)^ is equal to 77 (B) for 
all sets B such that (0, y) and (x, 0) doesn’t belong to B for any real x, y, 
that is, ^7^ (i), 7^ (t)^j represents the simultaneous jumps of the processes 
7j(t) and T 2 {t). 

Shortly speaking, we decompose the 2-dimensional CPP into the jump indepen¬ 
dent parts (T ^~and jump dependent parts (7^,7^). In other words, first 
part represents positive jumps only in one coordinate and second part repre¬ 
sents positive jumps in both coordinates. We denote by 11^,11^,1111 the Levy 
measures of the processes (7^,7^) resp., by A] 1 , A^, - the intensities 

of the corresponding processes, and by - the total number of jumps 

occuring only in the observed processes up to some fixed time T. 

The characterisitc function of two-dimensional CPP can be decomposed as 
follows (see [15]): 


E[exp(iziTi(t) + iz 2 T2(t))] = 

exp|t J { e M i^)-mi(d X ) +t J (exp(iz2y) — l)n^(dy)+ 

+t J (exp(iz\x + iz 2 y) — l)Illl(cte x dy) j- = 

E [exp(*zi7j‘ L (t))] E [exp(iz 2 T^(t))\ E exp(iziTi(t) + iz 2 T 2 (t)) 


In [17], Esmaeili and Kluppelberg applied the MLE approach to estimate the 
parameters in this model. According to the Theorem 4.1 from [17], the likelihood 
of the bivariate compound poisson process is given by: 

L(x 1 ,x 2 ,e 1 ,e 2 ,s) = h-i2-h 


(19) 
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where 

h = 

h = 

h = 


(Ai)” 1 exp(-A^T) 


n 

2=1 

(X 2 ) n * exp(-A^-T) 


fi{xi,Oi) (l „ F(u, \ 2 ', 6 )\ u=XiFi (~.. g ^ 


n 


h{yi,02) ( 1 q v F (Xi,v, 6 )\ v=x2 F 2 (y i -e 2 ) 


(AiA 2 ) n exp(—A n T) 


n 


d 2 


fi{xi-, 8 1 )h(yi-,e 2 ) ■ jr-^-F(u,v, 6 )\ u=Xl p l{x .. ei)> 

ouov v=\2F 2 ( yi -,e 2 ) 


where Xi and yt are jumps in the first and the second component occuring at 
different time, x and y are jumps for both components occuring at the same 
time, := / + °° fi(u;0i)du, i = 1,2. 

Formula (19) is based on a couple of simple facts, which directly follows from 
our construction. First, note that distance between jump moments for both 
components are independent and identically distributed exponential random 
variables with parameters Ai and A 2 . Second, one can show that 

Ui=:U± + ul ( 20 ) 

where Ui , * = 1,2- marginal tail integral, U^~ and uf are one-dimensional tail 
integrals of independent and dependent parts. Third, it is a worth mentioning 
that 


A 11 =C7(A 1 ,A 2 ;5). (21) 

From (20) and properties of two-dimensional tail integral and Levy copula we 
get that for any positive x, y 

\i(x)(l-Ff(x)) = ArA-C(AiA0r),A 2 ;«5)) 

A£(z)(l -i4'(a:)) = MF 2 {x)~ <7(Ai,A 2 F 2 (x);5)) 

A n F(x,y) = C(X 1 F 1 (x),X 2 F 2 (y)-6). 

For numerical example, we model the process of two-dimensional cumulative 
number of trades as a compond poisson process with exponential jumps. Denote 
by 0 \ and 8 2 the parameters of the jump sizes densities, that is, 

fi(x\ 6i) = 9i exp (—Oix) for x > 0, * = 1,2 

is the jump density for the i-th component. The dependence between 71 and 
71 is described by a Clayton Levy copula with parameter 5 . Clayton copula is 
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Table 1 


Estimated parameters 

Pair 

fh 

02 

5 

Ai 

A2 

Log likelihood value 

30-minutes returns 

Csco vs Int 

0,29 

0,14 

2,21 

24,91 

14,69 

5161,43 

Csco vs Msf 

0,23 

0,14 

2,71 

16,39 

17,21 

5196,93 

Int vs Msf 

0,14 

0,17 

2,38 

14,41 

24,68 

5579,32 

10-minutes returns 

Csco vs Int 

0,85 

0,43 

1,76 

74,18 

48,60 

10299,11 

Csco vs Msf 

0,71 

0,42 

2,11 

52,66 

55,78 

10406,96 

Int vs Msf 

0,42 

0,49 

2,00 

46,22 

72,48 

11511,90 


homogeneous copula and perfectly fits the conditions of the Theorem 5.1 (see 
Example 1 on page 7). 

The likelihood function of the continuously observed two-dimensional CPP 
process (71(f), 72 (f)) can be written in the following form, assuming that jumps 
occur at each moment for both components (there are no time intervals without 
trades in our data due to the fact that Cisco, Intel and Microsoft are liquid 
securities): 


L( Ai, A2, 0 \ 1 02 , < 5 ) 

( / n n 

= ((1 + 5)e 1 d 2 (\ 1 \ 2 ) s+1 ) n exp \ -AHr - (1 + <5) U + 02 


n (Ai exp(-0i<5a; l ) + A 2 exp(-0 2 %)) 


i =1 
— 1—2 


( 22 ) 


i= 1 


The results of the numerical optimization of In L(Ai, A 2 , 0i, 0 2 , S) are presented 
in the Table 1. 


6.3. Step 2. Estimation in time-changed model 

Let us shortly recall the model of stochastic time change in one-dimensional 
case. Denote by P t a equity price at moment t, and the returns by 


Y t = log 



(23) 


The main idea of the pioneer research [1] is to show that 


Yt = W(r(t)), 


(24) 


where W(t) = gt + oW(t) with Brownian motion Wt , and r(t) is the cumulative 
number of trades up to time t. It is a well- known fact that on one side, financial 
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Table 2 


Pair 

mi 

M2 



30-minutes data 

csco intc 

1.94E-08 

7,25E-09 

3,45E-09 

4,93E-09 

csco msft 

2,31E-08 

l,16E-08 

4,10E-09 

3,05E-09 

intc msft 

7,33E-09 

l,00E-08 

4,99E-09 

3,84E-10 

10-minutes data 

csco intc 

-6,00E-09 

-7,64E-10 

5.17E-10 

9,56E-10 

csco msft 

-7,05E-09 

-2,63E-09 

6,07E-10 

4,14E-10 

intc msft 

-7.85E-10 

-2,35E-09 

9,82E-10 

3,69E-10 


returns generally are not normally distributed (e.g. they are fat tailed), and on 
the other side, normality hypothesis is very convenient tool in finance, e.g. in 
mean-variance paradigm. In this context, formula (24) shows that returns are 
in some sense normal in business time, which differs from calendar time. 

Returning to our model (16)-(17), we now consider the problem of statistical 
estimation of the parameters p,\, ^ 2 , <Ji, er 2 . This task can separately solved for 
both components of the vector X(s) by the method of moments. Assuming that 
t is a CPP with intensity A and jumps distributed by exponential law with 
parameter 0 , we get 


E[Yt\ = M 


At 

T’ 


Var[F t ] 


a 2 \t t 2p, 2 \t 

~e ~ + e 2 


Solving the system of equations 

E[Y t } = E\Yt \, Var [F t ] = Va^], 


we arrive at the following estimates of the parameters ft and <r 2 : 

. _ 9E\Yt] ' 2 _ Var[At] - 2fi 2 \t/9 
^ At ’ ° At 

Estimation results are presented in the Table 1. 


6 . 4 ■ Step 3. Simulation techniques 

Here we show the performance of our approach introduced in Theorem 5.1. Our 
goal is to simulate two-dimensional time-changed Levy process: 

X(s) = (Wi(Ti(a)),W2(T2(s))) 

= (/xrTiOO + aiHMTiOO), M 2 T 2 (s)+u 2 ir 2 (T 2 (s))). 

Our simulation algorithm consists of the following steps: 

1. Model an independent sequence of jump times of a standard Poisson pro¬ 
cess r i: 

i 

Tj = ^ Tj < r, 

3 =1 
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where r determines the truncation level, Tj is a standard exponential ran¬ 
dom variable. 

2. Model k independent standard normal random variables and g[ 2 \ 
where i = 1 ,..., k. 

3. Model k independent uniform random variables Ri on [0,1], where i = 
1 ,--.,k. 

(3) 

4. Model k independent random variables G) ' with distribution function 
F(z), which is equal to 


F(z) 


(. z~ e + 1 ) 


(1 + 8 ) 


(25) 


by the method of inverse function, that is G,- 3 -* = F _1 (^), where are 
independent uniform random variables on [0,1]. 

5. Model subordinated Brownian motions by (truncated) series representa¬ 
tion: 


Z.is) := Y,^Ur 1 (T i )-G\ 1) -I{Ri<s}, 

i— 1 

Z 2 (s) := (V^Gf)) •/{#, < s}, 

i—1 * 

where the generalized inverse functions of Ui, i — 1,2 have the form 

UP'H*)={~k 1 ° g( J ) ’ f " l£A " (26) 

0, for x > A,;, 


and h*(x,y) is equal to h*{x,y) = xy. 

6 . Model two-dimensional subordinator (71 (s), 71 (s)) with Clayton-Levy cop¬ 
ula and compound poisson margins (with exponential jumps) by series 
representation (see [15], algorithm 6.13). 

7. Resulting trajectory is a linear transform of subordinator and subordi¬ 
nated brownian motion: 

*r(s) = thT^s) + alZ x {s) (27) 

X 2 (s) = fl 2 T 2 (s) + a 2 Z 2 (s) (28) 

Typical trajectories of simulated processes are presented in the Appendix. 
Figures 1 and 4 display trajectories for time-changed brownian motions mod¬ 
eled by Theorem 5.1 for 30 and 10 minutes data. Figures 2 and 5 show typical 
trajectories for suborninators modelled as compound poisson process with ex¬ 
ponential jumps for 30 and 10 minutes data. Finally, Figures 3 and 6 display 
resulting trajectories for the two-dimensional process X calculated by (27)-(28). 





V.Panov and I.Sirotkin/Series representations for time-changed Levy models 


18 


6 .5. Further research 

One interesting question, which was not addressed before, is to compare the Levy 
copula between simulated process ^Wi(71(s)), ^(^(s))^ and copula between 
subordinators (71 (s), 7-2 (s)). This question is motivated by the paper [16], where 
some relations between the corresponding correlation coefficients are given. 

In this paper, we would like to visually compare the copulas. The nonpara- 
metric estimation of the Levy copula between Levy processes X^ and X < ’ 2 ' > has 
been recently studied in [9]. The proposed estimator for any x\ > 0 and X 2 > 0 
is equal to 


F( Xl , x 2 ) = Y, 7{C/i,„(A£XW) < Xl , U 2 ,„(X n k X^) < z 2 }, (29) 

fc =1 

where by 

A n yW — x ^ — yW h — 1 n 7 — 1 9 

-^k A„ ~™~(k— 1)A„’ K 1..U, l 

we denote the increments of the processes X W and X^, and 

1 " 

U n 4x) = —Y 7{A£X (1) > x}, * = 1,2, (30) 


fc=i 


are the non-parametric estimators of the tail integrals of the underlined Levy 
processes. 

We applied this methodology to the simulated process ll'l (71 (s)) and W2(72(s)) 
and got the Levy-copula estimate (see Figure 7). 
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Appendix B: Graphs 




Fig 1. Time-changed brownian motion. Subordinators are CPP with exponential jumps. Pa¬ 
rameters are estimated from the Cisco and Intel 30-minutes data 




Fig 2. Subordinators for 30 minute data. Parameters are estimated from the Cisco and intel 
30-minute data 








first component & ~ first component first component 
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Fig 5. Subordinators for 10-minute data. Parameters are estimated from Cisco and Microsoft 
10-minute data 













